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Abstract 

The exchange coupling between Fe layers separated by BCC Cu is calcu- 
lated for Fe/Cu/Fe (001) trilayers. It is shown that the coupling is basically 
regulated by three extrema of the bulk BCC Cu Fermi surface. The con- 
tributions from those extrema are all of the same order of magnitude, but 
that associated with the "belly" at the T-point dominates. The calculated 
temperature dependence of the coupling varies considerably with spacer layer 
thickness. Individually, the amplitudes of these extrema contribution decrease 
with temperature, each according to a different rate. Such an effect may cause 

an actual increase of coupling with temperature for some Cu thicknesses. 
PACS:75.50.Fr., 75.30.Et., 75.50.Rr. 
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Although the common crystal phase of bulk Cu is FCC, it is possible to grow thin films 
of BCC Cu on Fe (001). The BCC stacking proceeds for up to 12 or 20 atomic planes 
approximately, but for larger thicknesses significant lattice modifications occur, leading to 
a structural transformation.!'! 

The exchange coupling between Fe layers separated by BCC Cu has been measured by 
groups at Simon Fraser University (SFU) and Philips.!'! Both have found that the coupling 
oscillates with decreasing amplitude as a function of the Cu thickness, but their results 
disagree in several important aspects.0 The Philips group data show well-defined short-period 
oscillations! whereas the SFU group originally observed a long-period oscillatory couplingJI 
Later, the SFU group found some indication of a short-wavelength oscillation in samples 
with smoother interfaces.0 The exchange coupling in multilayers can be strongly affected by 
sample interface quality! It is widely accepted that interface roughness tends to suppress 
short-wavelength oscillations and reduce the coupling amplitude. Therefore, as pointed out 
in Ref. [| it is rather puzzling that the values obtained at Philips are substantially smaller 
than those of the SFU group. 

Motivated by these apparently conflicting experimental results, we have undertaken a 
theoretical analysis of the exchange coupling between Fe layers across BCC Cu in Fe/Cu/Fe 
(001) trilayers. The coupling J is calculated for several temperatures T and spacer layer 
thicknesses N, using an extension of the formulation developed in Ref. [5[ For suficiently 
large N, we divide J(N) into oscillatory components coming from extrema which are re- 
lated to the spacer Fermi surface (FS). These oscillatory contributions to the coupling are 
calculated separately. Our results show that for perfectly smooth interfaces J(N) is domi- 
nated by short-period oscillations. We find that the temperature dependence of the coupling 
changes significantly with spacer layer thickness. The amplitude of each oscillatory compo- 
nent decreases with temperature, but they do so at different rates. We show that this may 
cause a surprising effect which is the increase of the coupling with temperature for some Cu 
thicknesses. 

The interlayer exchange coupling, defined as the total energy difference per surface atom 



between the antiferromagnetic and ferromagnetic configurations of the trilayer, is given byii: 

J = E/ duf(u)F(k h u,N) (1) 



where 

F = —Imtr ln[l + S\G c ^ irn — G^ m )S^((jJ +ljm+1 — Gm +1>?n+1 )], (2) 

S 1 ' = — G^yntG^+i m+itf)* 1 and = t(l — m +itfG'^ m t)~ 1 . In the equations above 
k\\ are the wave vectors parallel to the layers, f(u) is the usual Fermi-Dirac distribution 
function, and m is a plane index. As in Ref. |7] we consider an imaginary cleavage plane 
across the spacer between planes m and m + 1, separating the trilayer into two semi-infinite 
systems. G^ m and G^ +1 m+1 are matrices in orbital indices representing the surface one- 
electron Green functions of the left and right cleaved systems, respectively. The trace is 
taken over orbital indices and t denotes the spacer hopping matrix. 

This formula for the coupling is an extension of the result previously obtained in Ref. |5| 
and, for the one band model, reduces to the torque formula of Edwards et al.i In deriving it 
we have assumed that the electrons are noninteracting in the spacer and experience exchange- 
split one-electron potentials in the ferromagnetic layers. Most of the experimental results 
are for the bilinear exchange coupling term J\ which, for perfectly smooth Fe/Cu (001) 
interfaces, is virtually equal to J/2.! 

We have calculated the required Green functions within the tight-binding model with 
s, p, d orbitals and hopping to second nearest neighbours. The tight-binding parameters for 
all BCC Cu planes were determined from a first-principles LMTO-tight-binding electronic 
structure calculation of paramagnetic bulk BCC Cu. The parameters for ferromagnetic Fe 
were obtained from paramagnetic bulk Fei, by self-consistently adjusting the on-site energies, 
assuming charge neutrality. The effective intra-atomic electron-electron interactions were 
taken to be U^ p e = and U^ d = We neglect atomic potential differences due 

to the magnetic configuration change, thus making the approximation known as the "force 
theorem" . The kn sum in Eq.(l) is performed numerically and the energy integral is evaluated 



in the complex plane by summing over Matsubara frequencies at finite T. 

The calculated results at T=300K for the bilinear exchange coupling Ji as a function of 
Cu thickness are presented in Fig.|I] (full circles). Our results clearly show a short-period 
oscillatory exchange coupling, in excellent agreement with the Philips group data as far as 
the period of oscillation is concerned. 

For sufficiently large spacer thickness it is possible to express the coupling as a sum of 
oscillatory components whose periods are determined by extrema that are related to the 
spacer FS.EH1 It is essential to use a non-perturbative treatment, as in the quantum well 
approach, to analyze the relative importance of these contributions.0JlZl They depend upon 
the degree of confinement experienced by the carriers of both spin orientations in the corre- 
sponding extremum states, in the ferromagnetic and anti-ferromagnetic configurations of the 
magnetic layers.l'E10lii The widespread practice of considering only the periods predicted 
by RKKY theory, and treating the amplitudes and phases of these contributions as ad- 
justable parameters may be inadequate and rather misleading. The fitting usually involves 
several parameters and, in some cases, is not unique. Besides, when the spacer FS has to 
be regarded as consisting of more than one sheet, periods not predicted by RKKY theory 
may exist .0 Moreover, at finite temperatures, the decrease of the oscillation amplitudes 
as N increases is different for each extremum and may deviate strongly from the (1/iV 2 ) 
asymptotic regimei'0lll usually assumed in that sort of fitting. 

To identify the periods of oscillations of J(N) it is useful to look at the spacer FS. In 
BCC Cu only one energy band E(k\\, k^) crosses the Fermi energy Ep. Its calculated FS, 
shown in Fig.^(a), is basically a sphere with twelve "necks" developing at each face centre 
of the bulk BCC first Brillouin zone. In the (001) direction of growth, three sets of fcfj 
associated with the FS extrema contribute to the coupling. The first set consists of a single 
wave vector fcjj = (0, 0) (r-point) related to the FS "belly" . The other two are associated 
with the "necks". Set 2 consists of four vectors fcfj: (0, ±0.327) and (±0.327, 0), and set 3 of 
the M-points located at (±0.5, ±0.5). Here all wave vectors are given in units of 27r/a where 
a is the lattice constant. Due to the layered stucture of the system, it is useful to work with 
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the layer adapted bulk Brillouin zone instead of the usual BZ. The former is defined as a 
prism whose base is the two-dimensional first BZ and whose height is 2-K/d, where d is the 
interplane distance perpendicular to the layers. The relevant cross sections of the spacer FS, 
together with the corresponding extremal wave vectors /^(fcfj, Ep), are shown in Fig.Q(b). 
We must distinguish sets 1 and 2 from set 3 because, for the latter, the FS can be regarded 
as consisting of more than one sheet. This is because more than one extrema occurs in the 
first prismatic bulk BZ for each wave vector fcfj of set 3. 

Considering that the integrand F in Eq.|l| is an oscillatory function of N we can expand it 
in a Fourier series. However, it is necessary to generalize the expansion to a multiple Fourier 
serieslliil, when the equation E(k\\, Ar 1 ) = uj has more than one pair of solutions ±k^(k\\,uj). 
In this case, the general expansion of F is 

F(k h u,N)= £ ^.^LJ^**^™ (3) 

For N ^> 1 the exponential in Eq.||| oscillates rapidly as a function of k\\ and u. Thus, 
the stationary phase method can be applied, and the dominant contribution to the coupling 
comes from to = Ep and ku in the neighborhood of points at which the argument of the 
exponential is stationary. In this limit both the sum in k\\ and the energy integral in Eq.[l| 
can be evaluated analytically. The stationary points /cfj are the solutions of 

$>eVA^(%£ F ) = 0, (4) 

where V is the two-dimensional gradient in ku space. 

For fcjj belonging to sets 1 and 2 only one FS sheet occurs in the first prismatic bulk 
BZ; the analysis then proceeds exactly as in Ref. [7|. The corresponding periods are: p b =2.69 
atomic planes and p™=2.36 atomic planes respectively. However, for the M-points the FS can 
be regarded as consisting of two sheets. The two values of k x {k^ and k^) associated with 
/C|} M , shown by arrows in Fig. 0(b), correspond to equivalent periods p™ = 27r/2fc^ = 10.97 
atomic planes and p% = 27r/2fc^- = 1.1 atomic planes which cannot be distinguished just 
by looking at discrete values of N. However, Vkf- and V/c^ - vanish simultaneously when 
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calculated at (fc[j M , E F ). Thus, Eq. |] is satisfied for any values of n\ and n 2 , yielding other 
periods besides and p%- The relative contributions of these extrema depend on the com- 
parative values of the corresponding coefficients c nijn2 . The situation is very similar to that 
discussed in Ref. [15|. Nevertheless, our calculations have shown that the fundamental period 
(which is equivalent to p%) and its harmonics dominate. This is because the coefficients 
associated with them are far larger than those corresponding to alternative combinations of 
rii and n 2 . 

The calculated contributions to the coupling at T=300K coming separately from each 
set of fen are shown in Fig. [l|. We note that all three contributions are comparable, but the 
belly (full continuous line) clearly dominates. This contrasts with FCC Co/Cu (001) trilayers 
where the belly contribution is negligible in comparison with those given by the necks. 0The 
main reason for such difference is that minority carriers from the vicinity of the FS belly are 
fully confined in BCC Cu by the Fe layers in the ferromagnetic configuration. The physical 
origin of such confinement is that the sp-like BCC Cu band which intersects the FS has no 
counterpart at the minority-spin Fe FS. A similar situation happens for the Cu FS states 
in the vicinity of /cfj belonging to set 2, as shown in Fig.[| On the other hand, the Cu FS 
M-states of either spin can evolve into the corresponding Fe FS states because they have sp- 
character, due to the existence of a small but finite sp-d hybridization in this fc-space region. 
The degree of confinement experienced by the carriers in this case depends on the relative 
sp-d hybridization strengths. The agreement between the stationary phase approximation 
and numerical results verifies that the exchange coupling at room temperature in Fe/Cu 
(001) trilayers with perfect interfaces oscillates with a short period strongly influenced by 
the belly contribution. This is in accordance with the Philips group observations, as far 
as the period is concerned. However, the calculated strength is an order of magnitude 
larger than what they have observed. On the other hand the amplitude of our long period 
component is about three times smaller than the short period contribution. We believe that 
the discrepancy between experimental and theoretical results is due to interface roughness 
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which affects the amplitude and overall phase of the coupling.E2l The reason why the Philips 
results are much smaller than those of the SFU group remains unexplained. 

Motivated by recent measurements of the SFU group,0 we have calculated the temper- 
ature variation of the coupling for different Cu thicknesses. Our results, shown in Fig.[|, are 
for perfect interfaces where the bilinear exchange coupling J\ is much larger than the intrin- 
sic biquadratic term. The rate of variation of | Ji(T)\ with temperature clearly changes with 
spacer thickness. The calculated slope for Fe/12Cu/Fe agrees very well with the measured 
value for Fe/lOCu/Fe. The most striking result however, reproduced in the inset of Fig.|], 
is the increase of the coupling with increasing temperature for some Cu thicknesses. The 
temperature dependence of J\ is governed not only by the spacer FS but also by the con- 
fining strength of the ferromagnetic layersH, which differs for the three sets of extrema. As 



pointed out in Ref. [L9|, the energy dependence of the phases ip of the "Fourier" coefficients in 
Eq. |3| varies according to the confinement strength and is very important in determining the 
temperature dependence of the coupling. It turns out that the values of dip/duj calculated 
at the second set of extrema are about four and a half times larger than at the belly and the 
M-points. The temperature dependence of the former contribution is then stronger than the 
others. Hence, the coupling is approximately given by the sum of three oscillatory functions 
of N with comparable amplitudes which decay differently with temperature. Therefore, at 
some values of N, as in Fe/13Cu/Fe, the balance is such that an overall increase in the cou- 
pling is obtained even though the amplitude of each contribution separately decreases with 
temperature. Such increase was not detected by the SFU group. One possible explanation 
is that they have observed basically just a long period component. Another reason could 
be the influence of spin fluctuations in the ferromagnetic layers, which is neglected in our 
calculations. Nevertheless, this is an interesting temperature effect which may be observed 
under suitable conditions. 
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FIGURES 

FIG. 1. Calculated bilinear exchange coupling at T=300K for BCC Fe/Cu/Fe (001) trilayer 
versus Cu thickness (solid circles). The lines are contributions from the extremal points (see text) 
corresponding to: the Cu FS belly (full line), the neck wave vectors of set 2 (dashed line) and the 
neck M-points (dotted line). The inset shows the total contribution from the three sets of extrema; 
the tick are the same as those of the main figure. 

FIG. 2. Calculated BCC Cu FS (a), and its relevant cross sections for (001) (b). The arrows 
are the critical vectors Ar^fcfj). 

FIG. 3. Band structures of bulk BCC Cu and Fe in the relevant [001] direction for the wave 
vectors fc[j = (0,0) (a), fc[j = (0,0.324) (b) and fc[j = (0.5,0.5) (c). 

FIG. 4. Calculated temperature dependence of the bilinear exchange coupling for Fe/12Cu/Fe 
(solid circles) and Fe/14Cu/Fe (open circles). The inset is for Fe/13Cu/Fe; the tick labels are the 
same as those of the main figure. The lines are simply linear fits. 
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